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Abstract 

New method for finding exact solutions of nonlinear differential 
equations is presented. It is based on constructing the polygon corre- 
sponding to the equation studied. The algorithms of power geometry 
are used. The method is applied for finding one - parameter exact so- 
lutions of the generalized Korteveg - de Vries - Burgers equation, the 
generalized Kuramoto - Sivashinsky equation, and the fifth - order 
nonlinear evolution equation. All these nonlinear equations contain 
the term u m u x . New exact solitary waves are found. 

Keywords: the simplest equation method, travelling wave, exact solution, 
power geometry, the Korteveg - de Vries - Burgers equation, the Kuramoto 
- Sivashinsky equation, nonlinear differential equation. 

PACS: 02.30.Hq - Ordinary differential equations 

1 Introduction 



One of the most important problems of nonlinear models analysis is the 
construction of their partial solutions. Nowadays this problem is widely dis- 
cussed. We know that the inverse scattering transform [1-3] and the Hirota 
method [3-5] are very useful in looking for the solutions of exactly solvable 
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nonlinear equations, while most of the nonlinear differential equations de- 
scribing various processes in physics, biology, economics and other fields of 
science do not belong to the class of exactly solvable equations. 

Certain substitutions containing special functions are usually used for 
determination of partial solutions of not integrable equations. The most 
famous algorithms are the following: the singular manifold method [6-12], 
the Weierstrass function method [12-14], the tanh-function method [15-20], 
the Jacobian elliptic function method [21-23], the trigonometric function 
method [24,25]. 

Lately it was made an attempt to generalize most of these methods and 
as a result the simplest equation method appeared [26,27]. Two ideas lay in 
the basis of this method. The first one was to use an equation of lesser order 
with known general solution for finding exact solutions. The second one was 
to take into account possible movable singularities of the original equation. 
Virtually both ideas existed though not evidently in some methods suggested 
earlier. 

However the method introduced in [26, 27] has one essential disadvan- 
tage concerning with an indeterminacy of the simplest equation choice. This 
disadvantage considerably decreases the effectiveness of the method. In this 
paper we present a new method for finding exact solutions of nonlinear differ- 
ential equations, which greatly expands the method [26,27] and which is free 
from disadvantage mentioned above. When working out our method we used 
the ideas of power geometry recently developed in [28,29]. With a help of the 
power geometry we show that the search of the simplest equation becomes 
illustrative and effective. Also it is important to mention that the results 
obtained are sufficiently general and can be applied not only to finding exact 
solutions but also to constructing transformations for nonlinear differential 
equations. 

The paper outline is as follows. New method for finding exact solution 
of nonlinear differential equation is introduced in section 2. The application 
of our approach to the generalized Korteveg - de Vries - Burgers equation 
is presented in sections 3. Solitary waves of the generalized Kuramoto - 
Sivashinsky equation and of the nonlinear fifth-order evolution equation are 
found in sections 4 and 5, accordingly. 



2 Method applied 

Let us assume that we look for exact solutions of the following nonlinear 
n-order ODE 




(2.1) 
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In power geometry any differential equation is regarded as a sum of ordinary 
and differential monomials. Every monomial can be associated with a point 
on the plane according to the following rules 

C^yv — q 2 ), C 2 d k y/dz k — (-k, 1). (2.2) 

Here Ci,C 2 are arbitrary constants. When monomials are multiplied their 
coordinates are added. The set of points corresponding to all monomials of 
a differential equation forms its carrier. Having connected the points of the 
carrier into the convex figure we obtain a convex polygon called a polygon 
of a differential equation. Thus the nonlinear ODE (|2.1j) is associated with 
a polygon L\ on the plane. 

Now let us assume that a solution y(z) of the basic equation can be 
expressed through solutions Y(z) of another equation. The latter equation 
is called the simplest equation. Equations having known general solutions or 
solutions without movable critical points are usually taken as the simplest 
equations. Consequently we have a relation between y(z) and Y(z) 

y(z)=F(Y(z),Y z (z),...,z). (2.3) 

The main problem is to find the simplest equation. Substitution ()2.3)1 into 
the basic equation ({23} yields transformed differential equation which is 
characterized by new polygon L 2 . The suitable simplest equations should be 
looked for among the equations whose polygons, first, are of lesser or equal 
area than L 2 and, second, have all or certain part of edges parallel to those 
of L 2 . Let us suppose that we have found such polygon L 3 . Then we can 
write out the simplest equation 

E m (Y(z),Y z (z),...,z)=0. (2.4) 

It is important to mention that the choice of the simplest equation is not 
unique. If the following relation 

M n (F(Y, Y z ,..., z)) = RE m (Y, Y z ,...,z) (2.5) 

(where R is a differential operator) is true then it means that for any solution 
Y(z) of the simplest equation (|2.4J) there exists a solution of ()2.1|) . 
In this paper we will consider of identity substitution, i.e. 

y(z)=Y(z). (2.6) 

Then we should study the polygon L x of the equation ()2.1|) and try to 
find the polygon L3 generated by the simplest equation ()2.4j) . In this case, 
our method can be subdivided into four steps. 
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The first step. Construction of the polygon L±, which corresponds to the 
equation studied. 

The second step. Construction of the polygon L3 which characterizes the 
simplest equation. This polygon should posses qualities discussed above. 

The third step. Selection of the simplest equation with unknown param- 
eters in such a way that it generates the polygon L 3 . 

The fourth step. Determination of the unknown parameters in the sim- 
plest equation. 

To demonstrate our method application let us find exact solutions of 
nonlinear differential equations belonging to the following class 

n Ak I \ 

X^-CM + ^^W-O- (2-7) 

k=l 

This class of nonlinear ordinary differential equations containes a number 
of differential equations that corresponds to the famous nonlinear evolution 
equations. They are the Burgers - Huxsley equation, the Korteveg - de Vries 
- Burgers equation, the Kuramoto - Sivashinsky equation and so on. The 
following points M 1 = (-1,1), M 2 = (-2,1), . . ., M n = (-n,l), M n+1 = 
(0, 1), and M n+2 = (0, m+ 1) are assigned to the monomials of the equations 
(|2.7J) . The polygon L\ corresponding to the equations (|2.7J) is the triangle 
presented at figure 1 ((a): m > and (b): m < 0). 

Examining the triangle L\ we can easily find the polygon L3 assigned to 
the simplest equation of the first order (see smaller triangles at figure 1). 
Hence the simplest equation can be written as 

y z (z) = Ay^{z) +By(z). (2.8) 

This equation can be transformed to a linear equation. Its general solution 
takes the form 

, . / _ ( B mz} A\ m . 

y(z) = Ci exp — - - . (2.9) 



n J B r 

Later, having found the simplest equation, we should determine the co- 
efficients A and B for members of the class ()2.7|) . Then exact solutions of 
these equations will be expressed through the general solution of the simplest 
equation. 
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<?2 A 




(a) (b) 



Figure 1: Polygons corresponding to the equations studied and to the sim- 
plest equations. 

3 Exact solutions of the generalized Korteveg 
— de Vries — Burgers equation 

Let us find exact solutions of the generalized Korteveg - de Vries - Burgers 
equation with a help of our method. This equation can be written as 



u t - au m u x + (3u xxx = vu xx . (3.1) 

At m = 1 and v = (jH.lj) is the famous Korteveg - de Vries equation. 
Cauchy problem for this equation is solved by the inverse scattering transform 
[1]. Soliton solutions of this equation were found using the Hirota method [4]. 
At m = 1 and (/ / a dissipation is taken into account. In this case, the 
equation (jH.lj) is not integrable one. Its special solutions were obtained in 
the work [9]. Later these solutions were rediscovered many times. 
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Using the travelling wave reduction 

u(x, t) — y(z) , z = x — C t (3.2) 
and integrating with respect to z, we get 

(3y zz -vy z - C y — y m+1 = 0, (3.3) 

m + 1 

where the constant of integration is equated to zero. 

Taking into account the results obtained in section 2, we can write the 
simplest equation in the form 

y z (z)=Ay S! P(z)+By(z). (3.4) 
The general solution of the equation ()3.4j) is 

, W= ( Cl e X p{-^}-£)~*. (3.5) 

Let us formulate our results in the form of the following theorem. 
Theorem 3.1. Let y(z) be a solution of the equation 



, . v2 a m+2 , „ 2 v . . . 

Vz(z) = ± , V~*~(z) + -z-, sV(z) (3-6) 

K> y//3 (m + 2)(m + l) y VJ P(m + 4) yKJ 1 ; 

Then there exists a solution of the equation ()3.3j) that coincides with y(z) 
provided that m / 0, m / — 1, m ^ —2, m ^ —4, and 



2u 2 (m + 2) 
(3{m + Af 



Cc=-^™. (3.7) 



Proof. Let M(y 

i yzi Uzz) — be the generalized Korteveg — de Vries — Burgers 
equation ()3.3j) . By E(y,y z ) = denote the equation (|3.fij) . Substituting 

y z (z) = E{y,y z ) + Ay^(z) +By(z) (3.8) 

into the equation ()3.3j) and equating expressions at different powers of y(z) 
to zero we get algebraic equations for parameters A, B, and Cq in the form 



pA 2 (m + 2)(m + l)-2a = 0, (3.9) 
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[m 



+ l)(f3B 2 -uB - C ) = 0, (3.10) 



A (m + 1) ((3Bm + A(3B - 2 u) = 0. (3.11) 
Equations ()3.9|h ()3.10|) and (|3.11|) can be solved and we obtain 



A x 2 = ± , (3.12) 

v//?(m + 2)(m+l) 



C = -B (v-/3B), (3.13) 



B =jra- (3 - 14) 

Consequently we have found the parameters A, 5, Co of the theorem and 
the relation 

M{y,y z ,y Z2 ) = RE{y,y z ), (3.15) 
where R is a differential operator. □ 
Thus the solution of the equation ()3.3j) is found. It takes the form 



y(z) = — v M v ; + d exp <^ -— -N , 3.16 

1 ^2 z/ 2 (m + 2)(m + l) \ /? (m + 4) r 1 

where Ci is an arbitrary constant. 

Note that in the case, m — > the value of solitary wave velocity (J3.7)) 
tends to Co = — (4z/ 2 )/(16/3), but in the case, m->oowe get Co — > 0. 

It can be shown the correctness of the following formula 



(So + C, exp{-kz}y p = (2ff r p (l + tanh<; H^fA \]\ (3.17) 
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where So is an arbitrary constant and (fo is a constant connected with the 
coefficient C\ by the relation 



^ = iln^. (3.18) 
k S 

Using (|3.17|) the solution ()3.16|) can be rewritten as 



— — 
( v 1 (m + i)(m + 2)\ m ( {vm(z + ip a ) 1 \ m 

2/ = V T. A ; 1 + tanh <^ n ) ^° } . 3.19 

Assuming m — 1 in 1)3.19)1 . we get a solution of the Korteveg - de Vries - 
Burgers equation in the form 



This solution was found in [9]. 

4 Exact solutions of the generalized Kuramoto 
— Sivashinsky equation 

Let us look for exact solutions of the generalized Kuramoto - Sivashinsky 
equation, which can be written as 



u t + au m u x + 5u xx + I3u xxx + ju xxxx = 0. (4.1) 

Equation 1)4.1)) at m = 1 is the famous Kuramoto - Sivashinsky equation 
[30,31], which describes turbulent processes. Exact solutions of the equation 
1)4.1)1 at (3 = were first found in [30]. Solitary waves of 1)4.1)) at (3 ^ were 
obtained in [9] and periodical solutions of this equation at (3 ^ were first 
presented in [12]. 

Using the variables 



x 



x 



'6 , 5 2 
-, t' = t—, 



u 



u. 



a 



7 



7 



a 



5 2 



(4.2) 
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we get an equation in the form (the primes are omitted) 



u t + au m u x + u xx + auxxx + u xxxx = 0. (4.3) 

Taking into account the travelling wave reduction 
u(x, t) = y(z), z = x — Cot (4.4) 
and integrating with respect to z, we get the equation 

y zzz + <ry zz + y z - C y + — ^— y m+1 = 0. (4.5) 

m + 1 

Here a constant of integration is equated to zero. 
In this case, the simplest equation is the following 

m -j- 3 

y z = Ay— +By, (4.6) 

where A and B are parameters to be found. 
The general solution of ()4.6|) takes the form 

sW= ( Cl exppf£}-£p (4.7) 

Let us present our result in the following theorem. 
Theorem 4.1. Let y(z) be a solution of the equation 
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m+3 



Vz(z) — n n y 3 (z) t — f= y(z). 

1 2m 3 + llm 2 + 18m + 9; y V ; V2 m 2 + 18 m + 27 

(4.8) 



T/ten y(z) is a solution of the equation (j4.5j) in the case, m ^ ; m ^ —1, 
m ^ -§, m ^ -3, m ^ ~ 9± 2 3v ^ , and 

3(2m 2 + 9m + 9) 3 (3 + m) 

(2m 2 + 18m + 27) 3/2 v / 2m 2 + 18m + 27 



Proof. Let M(y,y z ,y zz ) = be the generalized Kuramoto - Sivashinsky 
equation (|4.5jl . By E(y,y z ) = denote the the equation Ij4.fij) . Substituting 

y*(z) = E{y,y z ) + Ay^(z) + By(z) (4.10) 
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into equation ()4.5|) and equating coefficients at powers of y(z) to zero, yields 
algebraic equations for parameters A, B, Cq and a in the form 



A 3 (2m + 3)(m + 3)(m + l) + 9a = 0, (4.11) 



+ 1) (B 3 + aB 2 + B- C ) = 0, (4.12) 



A 2 (3 + m)(m + l)(Bm + 3B + a) = 0, (4.13) 



(m + 1) (2(T 2 m 2 + 18a 2 m + 27a 2 - 9m 2 - 54m - 8l) = 0. (4.14) 
Solutions of the equations (|4~TTj) . (j4~T2|l . (|4~T3J) and (jQjj) can be written 

as 

A x = (-9 a) 1/s ( m 3 + 11 m 2 + 18 m + 9)~ 1/3 , 

/ (4-15) 
A 2j3 =(-9a) 1 / 3 f-i T iy (m 3 + 11 m 2 + 18 m + 9) " 1/3 , 



C = B (B 2 + aB + l), B=--^—, (4.16) 
v ' 3 + m 



a = ± - 3(3 + m) . (4.17) 
V2 m 2 + 18 m + 27 

In final expressions we will take into consideration only. 
So we have found parameters of the equation ()4.7|) . conditions ()4.9|) . and 
the relation 

M(y,y z ,y zz ) = RE{y,y z ), (4.18) 
where R is a differential operator. This completes the proof. □ 
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The general solution of the equation ()4.9|) can be presented in the form 



y (z) = [ C + Ct exp \±-= mZ -A) m , (4.19) 

where C± is an arbitrary constant and C is determined by expression 



c ± ^9^V2m2 + 18m + 27 
3{/(m + 3) (2m 2 + 5m + 3)' 

Again we note that the value of solitary wave velocity ()4.19)1 tends to 
Co = =F 1/V27 asm-> 0, but at the same time Co — > as m — > 00. 
Assuming m = 1 in (j4.19)l . we obtain the known solitary wave 



. V / 47v / -225a ^ (x-C t\\~ 3 „ 60 
2/0) = ^ + ^ ex P \ ± 7=^ f > = =F 



30 r L V47 J / 47V47 

(4.21) 



and the value of the parameter a: a = ±12/y47. This solution of the 
Kuramoto - Sivashinsky equation is the kink [9-13]. 

5 Exact solutions of the fifth — order nonlin- 
ear evolution equation 

Consider the fifth-order nonlinear evolution equation of the form 

Ut O. U U x -\- <J U xx + 5 U X xx ft u xxxx -\- u xxxxx 0. (5-1) 

It has the travelling wave solution 
u(x, t) — y(z) , z = x — Cot, (5.2) 
where y(z) satisfies 

Vzzzz + P Vzzz +Sy zz +cry z - C y — y m+1 = 0. (5.3) 

m + 1 
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Without loss of generality it can be set a — m+ 1. The simplest equation 
for ()5.3|) is the following 

y z (z) = Ay^iz) + By(z), (5.4) 

where A and B are parameters that should be found. 
The general solution of ()5.4j) is 

Let us summarize our results in the theorem. 
Theorem 5.1. Let y(z) be a solution of the equation 



yJz) = — = y * (z) — — — y(z). (5.6) 

y! ^/(m + 4)(3m + 4)(m + 2) KJ 3m + 8 yKJ V ; 



Then y(z) is also a solution of the equation (15. 3 j) provided that m ^ 0, 
|, m ^ -2, m ^ -§, 



m 7^ —1, m 7^ — ~ ; m ^ —2, m ^ — ~, m 7^ —4, and 



= /3 4 (m + 4)(3m + 4)(m + 2) 
° 2 (3m + 8) 4 1 ' ' 



P / 3 2 (llm 2 + 72m + 96) /3 3 (m 2 + 12 m + 16) , . 

= q , a = n . 5.8 

4(3m + 8) 2 4(3m + 8) 2 

Proof. Let M(y,y z ,y zz ) = be the equation (|5.3j) . By E(y,y z ) = denote 
the equation ()5.6|) . Substituting 



^(2) = + Ay^(z) + (5.9) 

into M(y,y z ,y zz ) = and equating expressions at different powers y(z) to 
zero, we obtain 



A 4 (m + 4)(3m + 4)(m + 2)-32 = 0, (5.10) 



C - B 4 - /3B 3 - 5B 2 - aB = 0, (5.11) 



12 



(m + 4)(m + 2)(3m J B + 2/? + 8S) = 0, (5.12) 



(m + 4) (7m 2 B 2 + 48mB 2 + 9QB 2 + 12(3mB + 48/35 + 16 5) =0, 

(5.13) 



— A m 3 B 3 + - A m 2 B 3 + — Am B 3 + 36 A B 3 + — 3 Am 2 B 2 + 
64 4 2 16 H 



(5.14) 



27 9 
+— 6 Am B 2 + 27 (3 A B 2 + -6 Am B + 18 5AB + 9aA = 0. 
4 4 



Solutions of the equations (I5~TUJ) . ffHTTTl) . (I5~T2"J) . (f5~T3j) . and ([534]) can be 

written in the form 



2 [2 

Ax 2 = ±— p (5.15) 
^2(m + 4)(3m + 4)(m + 2) 



A 3A = + 



(5.16) 



-2 v/2(m + 4) (3m + 4)(m + 2) 
C = (3B 3 + 5B 2 + aB + B\ (5.17) 



B = --?£-, (5.1* 
3m + 8' V 



5 = -X m 2 £ 2 - 3 mfi 2 -6B 2 -^0mB -30B, (5.19) 



I I 3 I 

a = m 3 B 3 - -m 2 B 3 - -mB 3 -4B 3 Bm 2 B 2 - 

64 4 2 16^ 



(5.20) 



-3mB 2 -3f3B 2 --5mB-25B. 
4 4 
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In final expressions Ai, A 2 will be used only. 

From obtained expressions (|5.17|) . f)5.18j) . (|5.19j) . and ()5.20|) we get con- 
ditions 1)5 .7|) . (|5.8j) . and the relation 

M{y,y z ,y zz ) = RE{y,y z ), (5.21) 
where i? is a differential operator. 

□ 

The general solution of f|5.6|) can be written as 



, , ^2(3m + 8) „ f /?mz 



^//5(m + 4)(3m + 4)(m + 2) ^i2(3m + 8) 

(5.22) 

where C\ is an arbitrary constant. 

Let us note that in the case, m — > the value of solitary wave velocity 
(522) tends to C = -/? 2 /256, but C -> as m -> cx). 



6 Conclusion 

In this paper a new method to look for exact solutions of nonlinear differ- 
ential equations is presented. Our goal was to express solutions of equation 
studied through general solutions of the simplest equations. The basic idea 
of our approach is to use geometrical representations for nonlinear differ- 
ential equations. The introduced method allows one to find the suitable 
simplest equations. Consequently our method is more powerful than many 
other methods that use an a priori expressions for y(z). In our terminol- 
ogy these methods take the a priori simplest equation. With a help of our 
approach we have found exact solutions of the following nonlinear differen- 
tial equations: the generalized Korteveg - de Vries - Burgers equation, the 
generalized Kuramoto - Sivashinsky equation and nonlinear fifth - order dif- 
ferential equation. All these equations belong to the class (12.7)) . Finally we 
would like to mention that this class can be also expanded. Exact solutions 
of equations belonging to the new class can be constructed with a help of our 
method. 
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